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Abstract 

An exact but simple general relativistic model for the gravitational 
field of active galactic nuclei is constructed, based on the superposition 
in Weyl coordinates of a black hole, a Chazy-Curzon disk and two rods, 
which represent matter jets. The influence of the rods on the matter 
properties of the disk and on its stability is examined. We find that 
in general they contribute to destabilize the disk. Also the oscillation 
frequencies for perturbed circular geodesies on the disk are computed, 
and some geodesic orbits for the superposed metric are numerically 
calculated. 

PACS numbers: 04.20.Jb, 04.40.-b, 98.58 Fd, 98.62 Mw 



1 Introduction 

There is a strong observational evidence that active galactic nuclei (AGN), 
X-ray transients and gamma-ray bursts (GRBs) are associated with accre- 
tion onto black holes, and that these sources are able to form coUimated, 
ultrarelativistic flows (relativistic jets). 

*e-mail: danielvt@ifi.unicamp.br 
^e-mail: letelier@ime.unicamp.br 
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The exact mechanisms to explain the production of jets are still uncer- 
tain, but they probably involve the interaction between a spinning black 
hole, the accretion disk and electromagnetic fields in strong gravitational 
fields (see, for example, and references therein). 

Thus, a reasonably accurate general relativistic model of an AGN would 
require an exact solution of Einstein-Maxwell field equations that describes a 
superposition of a Kerr black hole with a stationary disk and electromagnetic 
fields. Not even an exact solution of a stationary black hole-disk system has 
been found yet. Solutions for static thin disks without radial pressure were 
first studied by Bonnor and Sackfield J, and Morgan and Morgan and 
with radial pressure by Morgan and Morgan 0. Several classes of exact 
solutions of the Einstein field equations corresponding to static thin disks 
with or without radial pressure have been obtained by different authors 
[7]-jl6j. Thin rotating disks were considered in |171 118j . Perfect fluid disks 
with halos !19l and charged perfect fluid disks [20] were also studied, whereas 
thick relativistic disks were reported in j^J . Several solutions of the Einstein 
equations coupled to matter that represent disks has also been studied by 
the Jena group [221 "1201 • 

The static superposition of a disk and a black hole was first considered 
by Lemos and Letelier Zellerin and Semerak [201 found a stationary 

metric that reduces to the superposition of a disk and a black hole in the 
static limit and thus may represent a stationary disk-black hole system. The 
analysis of their solution is complicated by the fact that the metric functions 
cannot be analytically computed. For a survey on self gravitating relativistic 
thin disks, see for instance j31j . 

The aim of this paper is to consider the gravitational field of an AGN 
through a simple model: the static superposition of a black hole with a 
Chazy-Curzon disk and two rods placed on the symmetry axis, which will 
represent jets. Our principal interest here is to see how the presence of the 
rods affect the matter properties and stability of the disk. 

The article is divided as follows. In Sec. (21 we review the "displace, cut 
and reflect" method used to construct thin disks from a known solution of 
Einstein field equations in Weyl coordinates. Sec. [HI summarizes the for- 
malism to superpose thin disks and other Weyl solutions. Sec. [^ discusses 
Schwarzschild solution and the metric of a finite rod in Weyl coordinates. In 
Secinithe results of SeclSland^are then applied to construct the superposi- 
tion of disk, black hole and rods and the resulting energy-momentum tensor. 
In Sec. IHl the disk stability is studied through small horizontal and vertical 
oscillations about equatorial circular geodesies. In Sec. [7| some geodesic or- 
bits for the superposed metric are numerically calculated. Finally, Sec. IHlis 
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devoted to discussion of the results. We take units such that c = G = 1. 



2 Thin disk solutions in Weyl coordinates 

In absense of matter, the general metric for a static axially symmetric space- 
time in Weyl's canonical coordinates {t, r, z, (p) is given by 

= -e'^df + e^-<^(dr2 + dz^) + r'^e~'^dip^, (1) 

where (j) and v are functions of r and z only. Einstein vacuum field equations 
for the metric Eq. © yield [SI ES] 

r 

4',rT + — + <t>,zz = 0, (2a) 
r 

^[<t>\ = \ fr [{(t>l - 0'jdr + 2ct>,r(t>,zdz] . (2b) 



Given a solution of Eq. ()2a() - H2b|) . one can construct a thin disk by using the 
well known "displace, cut and reflect" method, due to Kuzmin First, a 
surface {z = 0) is chosen so that it divides the usual space in two parts: one 
with no singularities or sources, and the other with them. Then the part of 
the space with singularities or sources is disregarded. At last, the surface is 
used to make an inversion of the nonsingular part of the space. The result 
will be a space with a singularity that is a delta function with support on 
z = 0. The method is mathematically equivalent to make a transformation 
z ^ |2:| + a, where a is a constant. 

The application of the formalism of distributions in curved spacetimes to 
the Weyl metric Eq. is exposed in ,13^. One finds that the components 
of the distributional energy-momentum tensor [T^^] on the disk are 

-T\ = e'^-''{2-r4>^r)(f>,zS{z), (3) 
T% = e'^-^r<P^r<P,z5{z), (4) 
T% = T\ = 0, (5) 

where 5{z) is the Dirac distribution with support on the disk and is un- 
derstood that (p^z = lim2^o+ 4',z- The "true" energy density and azimuthal 
pressure are, respectively, 

a = e('^-'^)/2(-r^), (6) 
p = e^^-^^/^T^^. (7) 
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To explain the disk stability in absence of radial pressure, one may as- 
sume the counterrotating hypothesis, where the particles on the disk move in 
such a way that there are as many particles moving in clockwise as in coun- 
terckockwise direction. The velocity V of counterrotation of the particles in 
the disk is given by ES] 

V^ = ^ (8) 
a 

li < 1, the particles travel at sublumial velocities. The specific angu- 
lar momentum h of particles on the disk moving in circular orbits along 
geodesies reads 

h = T''l\-'^l\\- (9) 
V2(l-r</.,.) ' 

The stability of circular orbits on the disk plane can be determined with an 
extension of Rayleigh criteria of stability of a fluid at rest in a gravitational 
field: > 36 . Using Eq. © this is equivalent to 

<A,r(-3r</.,^ + 3 + r'^(t>^r) + r(t>,rr > 0. (10) 



3 Superposition of thin disks and other Weyl so- 
lutions 

An important property of the Weyl metric Eq. is that the field equation 
(|2aj) for the potential (p is the Laplace equation in cylindrical coordinates. 
Since Laplace's equation is linear, if and </>2 are solutions, then the su- 
perposition (/> = 01 + 02 is also a solution. The other metric function Eq. 
H2b|) is nonlinear, and so cannot be superposed. But one can show that the 
relation 

iy[ct>i + 4>2] = jy[4>i] + '^[^2] + 2iy[(l)i,(l)2], (ll) 

where 




- 4'l,z4>2,z)dr + {(I)l,r4>2,z + 01,^02,r)dz], (12) 



holds. Other useful relations are given in [S]. 

The energy-momentum tensor of the combined system disk and black 
hole has been computed by Lemos and Letelier ^HJ. Let cpn and (j)BH be 
the metric potentials of the disk and of the black hole, respectively. Then 
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the components [T^;,] of the superposition are 

-T\ = 6^^+^^"-" [2 - r{(t>D + 4>BH),r] (t>D,zS{z), (13) 

T'P^ = e't'^+'^^^-'^ri^D + (l>BH),r^D,z6{z), (14) 
T\ = T\ = 0, (15) 

where u = v[(f)£) + ubh], and again (t)D,z = li™2^o+ 4>D,z- The "true" energy 
density and azimuthal pressure read 

p= ^W-'PD-'f>BH)/2j^V^^ (17) 

Eq. (|13jl - (|14jl show that the potential of the black hole interacts with the 
disk and changes its matter properties. Although Eq. ()13() - H17|) have been 
derived for superposition of disk and black hole, they are also valid when 
the potential function (j)BH is a sum of other Weyl solutions, like the super- 
position of a black hole and rods. 



4 Black holes and rods in Weyl coordinates 

The Schwarzschild black hole metric function (pBH in Weyl coordinates is 
given by 

Vi+r2-2M\ , , 



')BH 



In 



ri+r2 + 2M, 

where r\ = (M — z)^ + and = (M + z)^ + r^. The function (j){r, z) can 
be related to the Newtonian potential U by 

(t) = 2U. (19) 

Thus, the metric potential of a finite rod of linear mass density A lying 
on the z axis and located along [ci,C2] is 



-2Aln 



C2 - z+ y^r2 + (c2 



ci - z + Y^r2 + (ci - z)^ 



(20) 



The calculation of the other metric function v for Eq. H18() . ()2U() and later 
for the superposed metric, is considerably simplified when one defines the 
following /i function 



IJ'k = Oik- z + \/r2 + {ak - zY, 



(21) 
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where is an arbitrary constant. This function is a natural consequence of 
the formaUsm of the inverse scattering method [371 138j . Eq. (jlBj) and ()2U|) 
can be rewritten as 



In ^ 
-2Aln 



(22) 
(23) 



where we defined 

= -M -z+ v^r2 + (M + z)2, 
//3 = C2 - z + \/r2 + (C2 - z)^, 

On using Eq. ((TT|) 



= M - z + + (M - z)2, 
Ci - z + \/r2 + (ci - z)2. 



/i4 



z^[ln /Xj — In /Xj] = i^[ln ^j] + z^[ln /ij] — 2z^[ln /ij, In ^j]; 



(24) 



the result 



i/[ln^j,ln/ij] = ln(/ij - /Xj), (25) 
which also follows from the inverse scattering method; and the identity 



(r^ + /xi/ij)(/xi - fJ-j) = 2(ai - aj)fiifj.j, 
one obtains following expressions for the metric function ly 

(r^ +;ui^2)^ 



i^BH = In 



i^i? = AX^ In 



L(^2+^2)(,2+^2)_ ' 
(r^ + ^3;U4)2 



(26) 

(27) 
(28) 



.(r2+^2)(^2+^2)_ 

5 Superposition of disk, black hole and rods 

We now consider the superposition illustrated in Fig. ^ a black hole with 
mass M whose center is on z = 0, two rods with equal mass density A, each 
one with mass M located along [— C2,— ci] and [ci,C2] on the z axis, and 
a disk on the plane z = constructed with the "displace, cut and reflect 
method" from the Chazy-Curzon solution with mass m, whose singularity 
lies on z = —a: 

0^ = (29) 

Y^r2 + {\z\ +a)2 
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M 



-M 
-<:2 

z=-a — • 

Figure 1: Superposition of a black hole with mass M, two rods and a Chazy- 
Curzon disk on the plane z = 0. 

It should be remembered that in Weyl coordinates a black hole with mass M 
is represented by a rod with length 2M, thus in Fig. ^we put a dotted circle 
around the rod in the middle. Such a configuration is not gravitationally 
stable: a consequence of the nonlinearity of Eq. ()2b|) is the appearence of 
gravitationally inert singular structures like struts between the rods and the 
black hole that keeps them appart; also in the superposition of the disk with 
the black hole, superlumial regions 

(y2 > 1) exist because there is matter 

up to the event horizon. 

The metric function cp of the superposition can be expressed as 

^ = _2Aln (^)+ln ( ^-^) - 2Aln ( ^) + cj^D, (30) 

with = -c\ - z + i/r2 + (ci + zY and /Ue = -C2 - z + y^r^ + (c2 + z)^. 

Now we consider the case when both rods just touch the horizon of the 
black hole, that is, when ci = M. Then = ^2 and ^5 = jii. From Eq. 



z=0 
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H13p - (|14() and Eq. (jH)), we get following conditions: 



cr > ^ y f2 + c| 

+ 2A(f2 + a2)3/2 
p>0- 



(f2 



\/l + > 0, 



f2 + c| - C2 



f2 + 



!yr^+(f2 +^2)3/2 



2A(f2 + a2)3/2 |^^^2 + g2 _ c2\/l + f2^ > Q, 



y2 < 1 



f2 + C2 



2otVi + ^^ + {f"^ + a^f''^{2 - \/l + f2) 



4A(f2 + a2)3/2/'^ 



' f2 + Co — C2 



\/l + f2^ < 0, 



(31) 



(32) 



(33) 



where f = r/M, a = a/M, £2 = C2/M, a = m/M, (3 = M/M and A = 
■M/{c2 — M) = /3/(c2 — 1). The conditions imposed are that of weak energy 
(a > 0), azimuthal pressure {p > 0) and sublumial velocity {V'^ < 1) of 
counterrotation of particles on the disk. For f — > 00, all three conditions are 
satisfied. In the regions where V"^ < 1, the weak energy condition is always 
satisfied, as can be seen by inequalities (|3T|) and (|33|) . 

Fig. [2Ia)-(c) shows curves of = 1 (dotted curves) and of hdh/dr = 
(solid curves) as function of A and f for three different values of a. The mass 
of each rod is kept constant (3 = 0.5 and the cut parameter takes values a = 1 
and a = 3 in (a) and (b), respectively. At the right of each dotted curve we 
have V'^ < 1 and the unstable regions of the disk appear between the curves 
of hdh/dr = 0. We note that in general with increasing mass of the disk and 
smaller length of the rods, the disk becomes more unstable and the regions 
of superlumial velocity also increase. There is, however, an interval of values 
for the rod's mass density where the zone of stability is increased, as can 
be seen in the lower part of the curve hdh/dr = for a = 0.001 in Fig. 
^a). This is probably due to the prolate quadrupole moment of the rods, 
which scale as AdP, where I is their length. Thus, for larger rods, the effect 
of prolate deformations may overwhelm the effect of the oblate quadrupole 
moment of the disk, and increase stability (see j39j for a detailed discussion 
of the effect of quadrupolar fields on the stability of circular orbits). 

Fig. IHta)-(b) shows again curves of = 1 (dotted curves) and of 
hdh/dr = (solid curves) as function of A and f for three different val- 
ues of a, but now the length of each rod is kept constant (62 = 11) and the 
cut parameter takes values a = 1 and a = 3 in (a) and (b), respectively. 
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2 4 6 8 10 12 14 16 




Figure 2: Curves ofV^ = I (dotted curves) and of hdh/dr = (solid curves) 
for the Chazy-Curzon disk in presence of a black hole and two rods. We keep 
the mass of each rod constant f3 = 0.5 and vary its mass density A. In (a)-(b) 
we take, respectively, a = 1 and a = 3. 
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Figure 3: Curves of = I (dotted curves) and of hdh/dr = (solid 
curves) for the Chazy-Curzon disk in presence of a black hole and two rods. 
We keep the length of each rod constant C2 = 11 and vary its mass density A 
(or equivalent, its mass). In (a)-(b) we take, respectively, a = 1 and a = 3. 




a=0.001 o 
a=0.5 • 
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With increasing masses of the disk and of the rods, the zones of instabil- 
ity and superlumial velocity of the disk are enhanced. In Fig. 0] (a)-(d) we 
plot the energy density a = Ma, azimuthal pressure p = Mp, square of 
the counterrotating velocity V and specific angular momentum h = Mh as 
functions of f for a = 3, a = 1, C2 = 11 (constant length) and different val- 
ues of the rod's linear mass density A. The curves were computed using Eq. 
(|13j) - (|17|) . Eq. 0-0 and Eq. (jSHI)- The expression for the corresponding 
metric function z/ is given in the Appendix. Energy density is lowered for 
a fixed radius as the rods become more massive, while pressure is slightly 
increased. Velocity of counterrotation and specific angular momentum are 
enhanced by increasing mass of the rods, as can also be deduced from Fig. 

Hb). 



6 Horizontal and vertical oscillations of the disk 

It is interesting to study the disk stability through the computation of hor- 
izontal (epicyclic) and vertical oscillation frequencies from perturbations of 
equatorial circular geodesies. Semerak and Zacek [301 have done such calcu- 
lations for the superposition of a Schwarzschild black hole with the Lemos- 
Letelier disk. They found that heavier disks are more stable with respect 
to horizontal perturbations near their inner rims, whereas they are less sta- 
ble with respect to vertical perturbations. For astrophysical relevance, it is 
important to determine not only the stability of circular motion on the disk 
plane, but also stability in the vertical direction. 

Using the perturbed equations for equatorial circular geodesies, the epicyclic 
frequency with respect to infinity iVh and the vertical oscillation frequency 
with respect to infinity iv^ for the metric are given by (see [101 ^ 
detailed deduction) 

p24,-u / 3 \ 

= n 7- ^,rr + rcj^l - 3<PI + -0,. , (34) 



2-r(h, 



[(t),,,-2(t>l{l-r(t>,r)\. (35) 



In Eq. ()35() the function (j)^ is obtained from the limit lim^^o* 4'z and 4>,zz 
follows from Eq. ((5aJ). Stable horizontal and vertical orbits are only possible 
where o;^ > and oj^ > 0, respectively. Note that condition > is 
equivalent to condition H1U() which follows from Rayleigh stability criteria. 

We compute first the frequencies for an isolated Chazy-Curzon disk, 
since is seems that such a calculation has not been done before for this 
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Figure 4: (a) Energy density a, (b) azimuthal pressure p, (c) eounter rotating 
velocity V'^ and (d) speeifie angular momentum h as functions of f for a = 1, 
o = 3, C2 = 11 and three different values for the rod's mass density A. 
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class of disks. In Fig. |S[a) we plot the epicyclic frequency as functions of 
radius r/m and cut parameter a/m. For a/m > 1.015 the disks always are 
stable and the epicyclic frequency is lowered for less relativistic disks. Highly 
relativistic disks (curve with a/m = 0.8 for example) develop annular regions 
of instability. The curves of vertical oscillation frequencies are depicted in 
Fig-Etb). We note that in this case highly relativistic disks are more stable 
in the vertical direction. In Eq. (|35|) the term with (pz is small compared to 
4)^zz , thus if we consider only 

2m .(^2_2„2)^ (3g^ 



(^2 + ^2)5/2 

we note that vertical oscillations are zero at r = a\/2, so the regions of 
vertical stability are enlarged as the cut parameter a is decreased. 

Now we consider the superposition of a Curzon disk and a black hole 
without rods. Fig. shows curves of (a) horizontal = Moj^ and (b) 
vertical uj^ = Muj^ oscillation frequencies of the disk with a = 1 and four 
different values of the "cut" parameter a. Now we always have regions of 
horizontal instability that begin at the innermost stable circular orbit and 
decrease for less relativistic disks. With respect to vertical oscillations, it is 
seen from Fig. El^b) that there are no regions of vertical instabilities. Thus 
one can conclude that the black hole desestabilizes the Curzon disk in the 
horizontal direction, whereas the opposite is true for the vertical direction. 

In Fig. [7|^a)-(b) we graph again iOh = Mwh and Coy = Mujy, respectively, 
for the Chazy-Curzon disk with black hole and rods, with their length fixed 
and vary the linear mass density A. As expected from the curves of Fig. 
I3fb), the more massive the rods, the larger are the disk's unstable regions in 
the horizontal direction. The rods also tend to lower the vertical oscillation 
frequencies near the disk's center, but unstable regions do not appear. 

7 Geodesic orbits 

In the previous section perturbations of equatorial circular geodesies were 
used to discuss disk stability for the system disk + black hole + rods. Now 
we solve numerically the geodesic equations of motion 

+ r^^x"x'3 = 0, (37) 

for metric Eq. where F^^ are the Christoffel symbols and the dot denote 
differentiation with respect to the proper time. Defining the orthonormal 
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Figure 5: Radial profiles of (a) the epicyclic frequency mujh and (b) the 
vertical oscillation frequency mujy for an isolated Chazy-Curzon disk. 
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tetrad {V'',W'',Y'', Z"} where 

= e"*/2(i, 0,0,0), (38a) 
PF" = e('^-'^)/2 (0,1,0,0), (38b) 
y'' = e(*-^)/2(o, 0,1,0), (38c) 

Z'' = (0,0,0,1), (38d) 

r 

the tetrad components of the four-velocity v"" read 

v"" = 7(l,wsin^/;cosX)^^sin'(/;sinx, t'cosV'), (39) 

with 7 = l/Vl — The specific energy and angular momentum of the test 
particle are 

£ = eH = e^^S, (40) 
h = r'^e'^ip = re'^^'^^v cos ip. (41) 

As initial conditions we take a position at radius ro on the disk's plane and 
components of the four- velocity Vq = 7(1, 0, sin^/;, cos V'), where vq is 
equal to the tangential velocity of circular orbits at radius tq. We choose 
initial radii such that the energy is slightly higher than the escape energy. 
Fig. IHta)-(b) shows the orbits of particles in the presence of the black hole 
and Curzon disk without rods. The parameters are a = 1, a = 3, fo = 3.9, 
£ ~ 1.01 and different initial angles ip. Fig. |HJa) is a projection of the orbits 
on the X — z plane. The coordinates have been transformed from Weyl to 
Schwarzschild coordinates {t, r, 6, ip) via the relations 

r = y^r(r - 2M)sin6l, z = (r - M) cos 6^, (42) 

and then to x = r sin 9 cos 99, y = r sin 9 sin (p and z = r cos 9. 

In Fig. inta)-(d) we have computed some orbits now with the rods. The 
parameters are q = 1, a = 3, fg = 7.43, A = 0.1, C2 = 11, £ ^ 1.01 
and different initial angles V- The orbit with ip = 89° has been placed in 
a separate graph for better visualization. For low initial angles, the rods 
have little effect on the trajectories, but this is not true as the particles 
approach the z axis. The orbit in Fig. intc)-(d) even suggests that we can 
expect chaotic behaviour for orbits that pass very near the rods. In fact, 
it has been shown 41 that prolate quadrupole deformations can introduce 
chaotic motion of geodesic test particles. In the oblate case, only regular 
motion was found. 
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Figure 8: Geodesic orbits for the superposition Curzon disk + black hole 
without rods. Parameters: a = 1, a = 3, f « 1.01, fo = 3.9. (a) Projection 
on the X — z plane of the curves in (b). 
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Ja) ^ ^ ^ , (b) 




Figure 9: Geodesic orbits for the superposition Curzon disk + black hole 
with rods. Parameters: a = 1, a = 3, A = 0.1, C2 = 11, £^ « 1.01 fo = 7.43. 
The curve for if) = 89° is displayed in (c) and (d) for better visualization. 
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8 Discussion 

We presented a very simplified, although exact, general relativistic model 
of an active galactic nuclei based on a superposition of a Schwarzschild 
black hole, a Chazy-Curzon disk and two rods placed on the symmetry axis, 
representing jets. We found that the presence of the rods enhances the disk 
regions with super lumial velocities. Using an extension of Rayleigh criteria 
of stability, it was found that in general the rods also increase the regions of 
instability, but when the rods are large and the disk's mass is low they can 
contribute to stabilize the disk. Also disk stability in the vertical direction 
was studied through perturbation of circular geodesies. The rods contribute 
to lower the vertical oscillation frequencies near the disk's center. Some 
geodesic orbits calculated numerically for the system black hole + disk + 
rods show the possibility of chaotic trajectories near the rods. 

The model here presented should be viewed as a first approach. As was 
stated in the introduction, more realistic models of active galactic nuclei 
should incorporate rotation and electromagnetic fields. However, the analy- 
sis of such a model would not be trivial, because of the large number of free 
parameters involved. 



D. V. thanks CAPES for financial support. P. S. L. thanks CNPq and 
FAPESP for financial support. 

9 Appendix 

The metric function Eq. H2b|) for the superposition of a black hole, two 
rods and a disk generated from Chazy-Curzon solution can be calculated as 
follows. We rewrite potential Eq. ()3U() as 
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4> = 4'Rl + (t>BH + 4>R2 + (p' 



(43) 



where 




with /i7 = — a 



e — z + Y^r^ + (—a — e — and fig 



= —a + € — z + 
bt. reduces to the 
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Chazy-Curzon disk Eq. ()29() on z > 0. Thus all terms can be expressed as 
fj, potentials. Using repeatedly properties H24() - H26() we get 

T^[(pRl + (/^BH + (pR2 + (Po] = + '^[(Pbh] + '^[(pR2] + 1^[</'J] 

+ 2l^[(l)Rl,(j)BH] +2iy[(j)Ri,4>R2] +2l^[(t)Ri,(p'jj] +2l^[(j)BH,4>R2] 
+ 2u[(t>BH,<Pi]+M<t>R2Ai], (44) 



with 



u[4>Ri] = AXHn 



i'[^R2] = 4A2 In 



(r^ + /i3^4)^ 

{r^+Hl^i2f 1 
_(r2 + /x2)(r2 + ^2)J ' 

(r2 +/i5^6)2 



(r2+^2)(^2+^2) 



2 2 

m r 



V 



J^[(pRi,(pBH] = 2 A In 

'Ri , 4'R2] = 4A^ In 

^[(l>Ri,<Pt)] = 7 7} , / , , ^ 

(a + ci)(a + C2)Vr2 + (a + 

-(a + ci)v^r2 + (c2 - z)^ + (ci - C2)^/r^+{a + z) 
J^[4>BH,4>R2] = 2Aln 



[r2 + (z + a)2]2' 

(r2 + ^i/X3)(r2 + ^2j»4) 
(r2 + ^i/X4)(r2 + ^2/U3) 

(r2 + /i3/x6)(r2 +^4^5) 
(r2 + /i3/X5)(r2 +/i4^6) 
2Am 



(a + C2)x/r2 + (ci 



(r2 + ^i^5)(r2 + /i2/»6) 

(r2 +^^^g)(r2 +^2/U5). 



-(a - M)v^r2 + (M - z)2 - 2MVr2 + (a + z)2j , 

i^[<Ai?,2, <Ad] = 7 77 ~T~7^^^=^f [(^2 - a)\/r2 + (ci + z)^ 

(a - ci)(a - C2)vr^ + (a + z)^ l 

-(ci - a)\/r2 + (c2 + 2;)2 + (ci - C2)\/r2 + (a + z) 
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In the particular case ci = M (//4 = /i2 and = //i), and on z = 0, Eq. 
H44() simplifies to 



'^[(f'Rl + (pBH + (t>R2 + 0j] = 
In 



+ 



C2(a2-M2)Jr2 + c2 



2 

2 



(r2 + a2)2 (a2 -c2)(a2 -M2)Vr2 + a2 
-M(a2 - c^)Vr2 + M2 + (M - C2)(a2 + Mc2)\/r2 + a 

(a2 - M2)Vr2 + a2 
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